Semi-classical analysis of real atomic spectra beyond Gutzwiller's approximation 
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Real atomic systems, like the hydrogen atom in a magnetic field or the helium atom, whose 
classical dynamics are chaotic, generally present both discrete and continuous symmetries. In this 
letter, we explain how these properties must be taken into account in order to obtain the proper (i.e. 
symmetry projected) h expansion of semiclassical expressions like the Gutzwiller trace formula. In 
the case of the hydrogen atom in a magnetic field, we shed light on the excellent agreement between 
present theory and exact quantum results. 

PACS numbers: 05.45.Mt, 03.65.Sq 



In the studies of the quantum properties of systems 
whose classical counterparts depict chaotic behavior, 
semi-classical formulas are essential links between the 
two worlds, emphasized by Gutzwiller's work [l|. More 
specifically, starting from Feynman's path formulation 
of quantum mechanics, he has been able to express the 
quantum density of states as a sum over all (isolated) 
periodic orbits of the classical dynamics. This formula, 
and extensions of it, have been widely used to understand 
and obtain properties of the energy levels of many clas- 
sically chaotic systems, among which the hydrogen atom 
in a maCTetic field |2, y{ , the helium atom |j, |^ |a| or 
billiardsl3,S,i,E3- 

At the same time, because the trace formula (and its 
variations) as derived by Gutzwiller only contained the 
leading term of the asymptotic expansion of the quantum 
level density, the systematic expansion of the semiclassi- 
cal propagator in powers of h has been the purpose of 
several studies [3, llOi ll^ ll^ i but which focused on bil- 
liards, for which both classical and quantum properties 
are easier to calculate. 

In a recent paper |l3|. general equations for efficient 
computation of h corrections in semi-classical formulas 
for a chaotic system with smooth dynamics were pre- 
sented, together with explicit calculations for the hydro- 
gen atom in a magnetic field. However, only the two- 
dimensional case was considered, because for the three- 
dimensional (3D) case, discrete symmetries and centrifu- 
gal terms had to be taken into account. Actually, this 
situation occurs in almost all real atomic systems depict- 
ing a chaotic behavior (molecules, two electron atoms...), 
for which experimental data involve levels having well 
defined parity, total angular momentum and, if relevant, 
exchange between particles. In particular, semi-classical 
estimations of experimental signals like photoionization 
cross-sections are calculated with closed orbits with van- 
ishing total angular momentum, whereas they usually 
involve P [L = 1) quantum states, whose positions in 
energy are shifted with respect to S* (L = 0) states. 
Furthermore, in recent years, the development of the 



harmonic inversion method makes it possible to extract 
the relevant quantities (position of peaks, complex am- 
plitudes) from both theoretical and experimental data 
with a much higher accuracy than with the conventional 
Fourier transform [lj|. In particular, it becomes possible 
to measure the deviation of the exact quantum results 
from the semi-classical leading order predictions. Thus, 
a detailed semi-classical analysis of experimental results, 
beyond the leading order in h, requires the understand- 
ing and the calculation of corrections due to both the 
discrete symmetries and centrifugal terms. In addition, 
we would like to stress that even if the present analysis 
is made with the density of states, it can also be made 
with the Quantum Green function, which leads to expres- 
sions and numerical computations of the first order h cor- 
rections for physical quantities like the photo-ionization 
cross-section [ij, |ig , which could either be compared to 
available experimental data [13, ll3l j or become a start- 
ing point for refined experimental tests of the quantum- 
classical correspondence in the chaotic regime. 

h corrections and discrete symmetries have already 
been discussed, but only for billiards ^ |l3, ll3l > whereas 
in the case of systems with smooth dynamics a detailed 
study is still lacking. Also, centrifugal terms and/or rota- 
tional symmetries have been considered by ma ny a uthors, 
but either in the case of integrable systems []j|,|2i|, or for 
values of the angular momentum comparable to the ac- 
tion of classical orbits ^1. 21 . 22] . From this point of view, 
the present study, which focuses on fixed values of the 
quantum angular momentum and the effect of the cen- 
trifugal terms on h corrections for systems with smooth 
chaotic dynamics, goes beyond the preceding considera- 
tions. More precisely, in this letter, we explain how to 
take into account both discrete symmetries and centrifu- 
gal terms in order to obtain a full semi-classical descrip- 
tion of the first order fi, corrections for the 3D hydrogen 
atom in a magnetic field. 

At first, in the case of a chaotic system, whose Hamil- 
tonian H = p^/2 + V{q) is invariant under a group S 
of discrete transformations a, the leading order of semi- 



classical approximation for the trace of the Green func- 
tion G{E) = 1/{E — H), restricted to the m*^ irreducible 
representation is given by [23j: 
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where the I sum is taken over all primitive (isolated) or- 
bits which become periodic through the symmetry oper- 
ation (7; (i.e. final position (resp. velocity) is mapped 
back to initial position (resp. velocity) by ai). Xmicrl^) 
is the character of aj'' in the rn"^ irreducible represen- 
tation of dimension dm- Si is the action of the orbit I, 
III the Maslov index, T/ the "period" , Af represents the 
Poincare surface-of-section map linearized around the or- 
bit and Ki is the subgroup of S leaving each point of the 
orbit I invariant. Adding first order h corrections, the 
preceding equation Q becomes : 



9rniE) = ^E^Ex™K)5S)(^)(l + ^<n) 



C*'jj can be derived by a detailed analysis of the sta- 
tionary phase approximations starting from the Feynman 
path integral, following the same steps as in Ref. 
and reads as follows : 
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leading to efficient computation of Gi,nit, t') and h correc- 
tions, that is, solutions of sets of first order differential 
equations, can easily be adapted to take into account 
these modified boundary conditions. 

As a numerical example, we have considered the 2D 
hydrogen atom in a magnetic field, at scaled energy 
e = —0.1 1^. More precisely, we have computed the trace 
of the quantum Green function, using roughly 8000 states 
belonging to the EEE representation |2J| of the group 
D4, corresponding to effective 1/fi values ranging from 
to 124 (See Ref. [13 for further details). In that case, the 
periodic orbit T234 ^^ (see inset of the top of Fig. [J 
for the trajectory in semi-parabolic coordinates), being 
(globally) invariant under a rotation of angle 7r/2, gives 
rise to contributions in the semi-classical approximation 
of the trace at all multiples of Sj234:/^- In the same way, 
the periodic orbit 1243 (see middle inset of Fig. ^ being 
invariant under a rotation of angle tt, contributions are 
present at all multiples of 81243/'^- For both these orbits, 
tabled displays the comparison of the present theoretical 
calculation and the numerical coefficient Cf^^ , extracted 
from the exa ct q uantum Green function, using harmonic 
inversion [l3lll4l |. As one can notice, the agreement is ex- 
cellent for the amplitudes and rather good for the phases, 
which is the usual behavior of harmonic inversion. Fur- 
thermore, the same agreement has also been found for 
the other representations, thus emphasizing the present 
approach for the calculation of the first order h correc- 
tions when taking into account discrete symmetries. 
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where Cf^ 



*^ arises from the time to energy domain 
transformation. Ci^n{to) (see Ref. 131 ^^ the expres- 
sions) involves the classical Green functions Gi^n{t,t'), 
i.e. the solutions of the equations controlling the linear 
stability around the classical trajectory q^'„(0 : 
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The fact that the orbits are periodic after the symmetry 
transformation af" determines the boundary conditions 
that the classical Green functions Qi^n{t,t') must fulfill, 
namely : 



TABLE I: Numerical comparison between the theoretical h 
corrections C\^ for the trace of the quantum Green function, 
restricted to the EEE representation, of the 2D hydrogen 
atom in a magnetic field and the numerical coeflicients C/^^ 
extracted from exact quantum function using harmonic in- 
version. The agreement is excellent for the amplitudes and 
rather good on the phases, thus emphasizing the validity of 
the present approach. 
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cr 


\pHI\ 


Rel. error 


argC«^ 


il234 
|l234 


-0.094 430 


0.09445 


« 2 X 10"" 


0.9996x7r 


-0.361689 


0.3611 


Ri 2 X 10"^ 


0.996x7r 


fl234 


-0.400 555 


0.3992 


Ri 3 X 10"^ 


1.005x7r 


il243 


0.049 339 


0.0493 


f« 8 X 10-^ 


-0.075x7r 



[ Qt^a^'QiAnTiX) = Qt„gi,„(0,t') 

(6) 
where Vto is the projector along the "periodic" orbit at 
the position depicted by time io and Qt„ = 1 — Vta ■ Of 
course, for cr; = 1, one recovers the boundary conditions 
given in Ref. Il3l|. Finally, all technical steps of Ref. |l3| 



Contrary to the preceding, calculating first order h cor- 
rections due to centrifugal terms is nrore complicated and 
is best explained in the case of the 3D hydrogen atom in 
a magnetic field. The regularized Hamiltonian in semi- 
parabolic coordinates, for fixed value M of the projection 
of the angular momentum along the field axis, is given 



H = 
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Hq is then the Hamiltonian of the 2D hydrogen atom in a 
magnetic field. If U{u, v) was regular, then the additional 
first order H correction for the orbit I would simply be : 

rT, 

dtU{ui{t),Vi{t)). (8) 
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transformation 1271 



4/^0 

that in this case, the Langer 

of the coordinates {u, v) — > 

(exp (—cc), exp (—?/)) gives rise to a Hamiltonian which 
does not separate into kinetic and potential energies and 
for which no expressions for h corrections are available. 
On the other hand, the fact that U{u, v) is singular 
imposes boundary conditions on both classical and quan- 
tum dynamics. The classical trajectories have to make 
(smooth) bounces near u = and v = and for van- 
ishing values of h, we expect the trajectories of H to be 
those of Hq, but mapped onto the reduced phase space 
{u > 0,v > 0), i.e. making hard bounces on the {u,v) 
axis. From the quantum point of view, depending on 
the parity of M, only wavefunctions belonging to given 
representations of Z?4 are allowed. Thus, first order h cor- 
rections due to the singular part of the potential U, are 
given by the preceding considerations on the symmetries, 
whereas remaining corrections are given by Eq. (jHJ , where 
U has to be replaced by a smooth counterpart, namely : 
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(9) 
Actually, one can show that the preceding equation gives 
the right answers for h expansion of the propagator of 
the free particle (up to h^) and the harmonic oscillator 
(up to ft^), for which analytical expressions for classical 
trajectories, classical Green functions and quantum prop- 
agators exist (higher orders have not been checked yet). 
However, even if a detailed analysis of the derivation of 
the trace formula in presence of centrifugal terms seems 
to show that the preceding approach works in general 
cases, rigorous proof of Eq. ^ is lacking. 

Nevertheless, in the case of the 3D hydrogen atom in a 
magnetic field, we have compared the first order h correc- 
tions, for different periodic orbits and for different values 
of the magnetic number M , with the present prediction, 
namely : 



Cr(Af) = Cr(2i?)--(4|M|2-l) 



dtU{ui{t),vi{t)). 
(10) 
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FIG. 1: First order h correction to the semi-classical approxi- 
mation of the trace of the quantum Green function for the hy- 
drogen atom in a magnetic field for different values of the rnag- 
netic number M, M — 1/2 corresponding to the 2D case p^ . 
Crosses depict the values extracted from the exact quantum 
function using harmonic inversion, whereas the solid line cor- 
responds to the classical results given by Eq. HlOjI . For the 
three different periodic orbits, whose trajectories in the {u, v) 
plane are plotted (the nucleus being depicted by the black 
dot), the agreement is excellent, thus emphasizing the valid- 
ity of Eq. @ and Eq. itTO . 



The results are displayed in Fig. Q] for M — 0, 1,2 and 
for three different orbits, namely 1234, 1243 and 12343, 
whose trajectories in the (u, v) plane are plotted. The 
solid line is the theoretical result given by Eq. H10|l . 
whereas the crosses are the values extracted from the 
trace of the exact quantum Green function, using har- 
monic inversion (for scaled energy e = —0.1, roughly 8000 
effective 1/h values ranging from to 124). As one can 
notice the agreement is excellent, thus giving strong sup- 
port for the validity of Eq. ^ and Eq. I|10() . Furthermore, 
the simplicity of the replacement U may serve as a guide- 
line for a rigorous treatment of the h corrections arising 
from the centrifugal terms. In particular, the calcula- 
tion of higher orders involves products of the derivatives 
of these centrifugal terms and those of the potential Vq, 
giving rise to non-trivial mixing between centrifugal and 
standard h corrections. 

In conclusion, we have presented a semi-classical anal- 



ysis, beyond the usual Gutzwiller approximation, includ- 
ing first order h corrections, of the quantum properties 
of real chaotic systems. More specifically, we have ex- 
plained the additional corrections arising when taking 
into account both discrete symmetries and centrifugal 
terms. In the case of the (3D) hydrogen in a magnetic 
field, the agreement between the theory and the numer- 
ical data extracted from exact quantum results is excel- 
lent, emphasizing the validity of the analysis, especially 
of equations © and H1Q|I . 

Finally, since we know how to compute the h correc- 
tions, it would be very interesting to work the other way 
round, that is, to perform the semi-classical quantiza- 
tion, thus getting h corrections in the semi-classical es- 
timations of the quantum quantities, like the eigenen- 
ergies. Of course, this represents a more considerable 
amount of work, since the C*'„ coefficients must be com- 
puted for all relevant orbits and then included in stan- 
dard semi-classical quantization schemes, like the cycle 
expansion [iHTilEil. 
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